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Abstract The huge datasets collected at the two B factories, Belle and BaBar, have made 
it possible to explore the radiative penguin process b — >• sj, the electroweak penguin process 
b — s> s£ + £~ and the suppressed radiative process b — >• dj in detail, all in exclusive channels and 
inclusive measurements. Theoretical tools have also advanced to meet or surpass the experi- 

p^ mental precision, especially in inclusive calculations and the various ratios of exclusive channels. 

^vi , In this article, we review the theoretical and experimental progress over the past decade in the 

radiative and electroweak penguin decays of B mesons. 
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I INTRODUCTION 

The B meson system, which is a bound state that consists of a b quark and a 
light antiquark, provides an ideal laboratory for precise study of the Standard 
Model (SM) of particle physics, and thus facilitates the search for new physics 
(NP). Because the b quark mass is much larger than the typical scale of the 
strong interaction, the otherwise troublesome long-distance strong interactions 
are generally less important and are under better control than in other lighter 
meson systems. Radiative penguin_| decays of the B meson with the emission of 
a photon (7) and electroweak penguin decays with the emission of a lepton pair 
(l + £~ , £ = e, /i) are of particular interest in this respect. These processes reveal 
the flavor changing neutral current (FCNC), that is the transition of a b quark 
with an electric charge of —1/3 into an s or a d quark of the same charge. In 
the SM, the FCNC occurs only via virtual loop diagrams (Fig. [[]). Additional 
NP contributions to these decay rates are not necessarily suppressed with respect 
to the SM contribution. Examples of such NP scenarios include those in which 
the SM particles in the loop diagram are replaced by hypothetical new particles 
at a high mass scale; so far, they have not been directly accessible in collider 
experiments. Radiative and electroweak penguin decays are highly sensitive to 
NP because they are theoretically well-understood and have been extensively 
measured at the B factories. The search for such NP effects complements the 
search for new particles produced at collider experiments. 

The first generation of the B factories at KEK (the Belle experiment at the 
KEKB e + e" collider) © and at SLAC (the BaBar experiment at the PEP- 

II e + e~ collider) (j3]) have collected huge samples of B meson decays and have 
thereby established the SM picture of CP violation and other flavor-changing pro- 
cesses in the quark sector. These processes are governed by a single 3x3 unitarity 
matrix referred to as the Cabibbo-Kobayashi-Maskawa (CKM) matrix (5,6). The 
CKM matrix can be illustrated by a unitarity triangle in the complex plane that 
is overconstrained by measurements from the B factories, the Tevatron B physics 
programs (namely the CDF ([7]) and DO (jHJ) experiments), and earlier kaon de- 
cay experiments. In other words, none of the current measurements of B meson 
decays have observed any unambiguous sign of NP (I MTfJj) . Although this exper- 
imental result is an impressive success of the CKM theory within the SM, there 
is still room for sizable new effects from new flavor structures, given that FCNC 
processes have been tested up to only the 10% level. 

The nonexistence of large NP effects in the FCNC processes hints at the famous 
flavor problem, namely why FCNCs are suppressed. This problem must be solved 
in any viable NP model. Either the mass scale of the new degrees of freedom is 
very high or the new flavor-violating couplings are small for reasons that remain 
to be found. For example, assuming generic new flavor-violating couplings, the 



lr The name penguin decays was first introduced in Ref. |T} as the result of a bet. A more 
detailed account of the name can be found in Ref. ([2]). 
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present data on K-K mixing implies a very high NP scale of order 10 3 -10 TeV 
depending on whether the new contributions enter at loop-level or at tree-level. 
In contrast, theoretical considerations on the Higgs sector, which is responsible 
for the mass generation of the fundamental particles in the SM, call for NP at 
order 1 TeV. As a consequence, any NP below the 1-TeV scale must have a 
nongeneric flavor structure. The present measurements of B decays, especially of 
FCNC processes, already significantly restrict the parameter space of NP models. 
For further considerations on NP, the reader is referred to another article in this 
volume (HH and to Ref. (fT2l) . 

Quark-level FCNC processes such as b — > sj, to which NP may contribute, 
cannot be directly measured because the strong interaction forms hadrons from 
the underlying quarks. Instead, the experimentally measured and theoretically 
calculated process is a B meson decay into a photon plus an inclusive hadronic 
final state X s , which includes all the hadron combinations that carry the strange 
quantum number s = +1 of the s quarko Exclusive final states with one or a 
few specific hadrons in the final state (e.g., B — > K*"f) have less predictive power 
theoretically; however, because the measurements are easier and better defined, 
there are other useful observables beyond branching fractions, in particular CP, 
forward-backward, isospin, and polarization asymmetries. In the future, a large 
overconstrained set of measurements of these observables will allow us to detect 
specific patterns and to distinguish between various NP scenarios. 

This review covers progress in radiative and electroweak decays in the past 
decade, during which a huge number of B factory results were accumulated and 
significant progress in various theoretical aspects was achieved. The pioneering 
work that led to the first observation of the 6—7-57 process by CLEO fj!3|.ll4j) 
was discussed in an earlier volume of this journal (|15p . 

Our review is organized as follows. In Section [2j we describe theoretical tools 
for radiative and electroweak penguin decays, and in Section[3]we describe exper- 
imental techniques. We give theoretical predictions in Section [4] and summarize 
the measurements of radiative and electroweak penguin decays in Section [5j Fi- 
nally, we briefly discuss future prospects in Section [H 

2 THEORETICAL FRAMEWORK 

Inclusive B decays are theoretically clean because they are dominated by par- 
tonic (perturbatively calculable) contributions. Nonperturbative corrections are 
in general rather small (|16 1fT71 18). This result can be derived with the help of 
the heavy mass expansion (HME) of the inclusive decay rates in inverse powers 
of the b quark mass. Up-to-date predictions of exclusive B decays are based on 
the quantum chromodynamics (QCD)-improved factorization (QCDF) and soft 
collinear effective theory (SCET) methods. In general, exclusive modes have 
larger nonperturbative QCD corrections than do inclusive modes. 



In this review, we use the following notations and conventions: We denote the inclusive 
decay as B — ¥ X s j when charge conjugation is implied, or as B — ► X s j and B — >• X-g-y to reflect 
the quark charges of the underlying processes b — ► sy and b —$■ sj, respectively, when CP and 
angular asymmetries are concerned. Here, B denotes either an isospin- and CP-averaged state 
of B , B°, B + and B~ mesons, or an isospin averaged state of B° and B + (in the latter case, 
B denotes B° and B~). Expressions are constructed similarly for inclusive XdJ and X s £ + £~ 
final states, and isospin-averaged exclusive final states. In the literature, the notation B — s* X s j 
is also commonly used for the case that includes charge conjugation. 
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2.1 Electroweak Effective Hamiltonian 

Rare B decays are governed by an interplay between the weak and strong inter- 
actions. The QCD corrections that arise from hard gluon exchange bring in large 
logarithms of the form aj(mj) log" 1 (m^/M), where M = mt or M = mw and 
m < n (with n = 0, 1,2,...). These large logarithms are a natural feature in any 
process in which two different mass scales are present. To obtain a reasonable 
result, one must resum at least the leading-log (LL) series, n = m, with the 
help of renormalization-group techniques. Working to next-to-leading-log (NLL) 
or next-to-next-to-leading-log (NNLL) precision means that one resums all the 
terms with n = m + l or n = m + 2. A suitable framework in which to achieve 
the necessary resummations of the large logarithms is an effective low-energy 
theory with five quarks; this framework is obtained by integrating out the heavy 
particles, which in the SM are the electroweak bosons and the top quark. 

This effective field theory approach serves as a theoretical framework for both 
inclusive and exclusive modes. The standard method of the operator product 
expansion (OPE) (19,20) allows for a separation of the B meson decay amplitude 
into two distinct parts: the long-distance contributions contained in the operator 
matrix elements and the short-distance physics described by the so-called Wilson 
coefficients. 

The electroweak effective Hamiltonian (21,22,23) can be written as 



n eS = ^-J2 c ^ M ) °*(m), (i) 

where Oi{u) are the relevant operators and Ci(//, M) are the corresponding Wil- 
son coefficients. As the heavy fields are integrated out, the complete top and 
W mass dependence is contained in the Wilson coefficients. Clearly, only within 
the observable H e s does the scale dependence cancel out. Gf denotes the Fermi 
coupling constant. 

From the [/, independence of the effective Hamiltonian, one can derive a renor- 
malization group equation (RGE) for the Wilson coefficients Cj(//): 

V-j-C i (u)='y ji C j (n.), (2) 

where the matrix 7 is the anomalous dimension matrix of the operators Oi, which 
describes the anomalous scaling of the operators with respect to the scaling at 
the classical level. At leading order, the solution is given by 



Qin) 



a s (nw) 



OL s {u) 



2/3 _ 



l + /3o^hi^ 



7?. 

2/3 

Ci(/J,w), (3) 



where \i d/d^i a s = — 2/3oa^/(47r), and j3q and 7° correspond to the leading anoma- 
lous dimensions of the coupling constant and of the operators, respectively. The 
tilde indicates that the diagonalized anomalous dimension matrix is used. 

Although the Wilson coefficients Cj(/x) enter both inclusive and exclusive pro- 
cesses and can be calculated with perturbative methods, the calculational ap- 
proaches to the matrix elements of the operators differ between the two cases. 
Within inclusive modes, one can use the quark-hadron duality to derive a well- 
defined HME of the decay rates in powers of A/m b d2H [251 EH [271 EH [29D . In 
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particular, the decay width of the B — > X s j is well approximated by the partonic 
decay rate, which can be calculated in renormalization-group-improved pertur- 
bation theory (|30ll3T]) : 



T(B ->■ X sl ) = T(6 -► X s parton 7) + A nonpert - . (4) 

In exclusive processes, however, one cannot rely on quark-hadron duality, so one 
must estimate the matrix elements between meson states. A promising approach 
is the QCDF-method, which has been systematized for nonleptonic decays in 
the heavy quark limit ([32J I33 J I34" )) . In addition, a more general quantum field 
theoretical framework for QCDF, known as SCET, has been proposed (135 ( 136 1 137 } 
38,39,40). This method allows for a perturbative calculation of QCD corrections 
to naive factorization and is the basis for the up-to-date predictions for exclusive 
rare B decays. However, within this approach, a general quantitative method to 
estimate the important A/nib corrections to the heavy quark limit is missing. 

2.2 Perturbative Corrections to Inclusive Decays 

Within inclusive B decay modes, short-distance QCD effects are very important. 
For example, in the B — > X s ^f decay these effects lead to a rate enhancement by 
a factor of greater than two. Such effects are induced by hard-gluon exchanges 
between the quark lines of the one-loop electroweak diagrams. The corresponding 
large logarithms have to be summed as discussed above. 

The effective electroweak Hamiltonian that is relevant to b — > s/d'y and b — > 
s/d(, + l~ transitions in the SM reads 



V - 4Gf 



10 2 



(5) 



where the explicit CKM factors are A^ = V tb V t * and A^ = V ub V*. The unitarity 
relations A^ = — A^ — A^ have already been used. The dimension-six operators 
are 

01 = (s L ^T a c L )(c L ^T% L ) , 2 = (sLl„c L )(c Ll »b L ) , 

0\ = (s Ll ^T a u L )(u L ^T a b L ) , Ol = (s Llfl u L )(u L ^b L ) , 

03 = {SLl^ L )Y, q {q^q) , 04 = (sL^T% L )j: q (qrT a q) 



05 = (s L Tb L )j: q (qr'q), 6 = (s L TT a b L )j: q (qT'T a q) , 

°7 = ji^rn b (s L a^b R )F^ , O s = ^ r m b (8 L a' tv T a b R )G^, , 

09 = it(^7A)E^7^) ; 0io = I f^r(sL7A)£^7 / W), 



(6) 



where T = j^ju'yx and F' = 7^7"7 . The subscripts L and R refer to left- and 
right-handed components, respectively, of the fermion fields. In b — > s transitions 
the contributions proportional to A" are rather small, whereas in b — > d decays, 
where AJ| is of the same order as A^; these contributions play an important role 
in CP and isospin asymmetries. The semileptonic operators 0$ and 0\$ occur 
only in the semileptonic b —> s/d£ + l~ modes. 

Among the four-quark operators, only the effective couplings for i = 1,2 are 
large at the low scale [i = mj, [Cxiim^) ~ 1], whereas the couplings of the other 
four-quark operators have almost negligible values. But the dipole operators 
[C-j{m b ) ~ —0.3, Csim^) ~ —0.15] and the semileptonic operators \C%{m b 
4, Cio(mt,) ~ —4] also play a significant role. 



r^J 
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There are three principal calculational steps that lead to the LL (NNLL) result 
within the effective field theory approach: 

1. The full SM theory must be matched with the effective theory at the scale 
H = fxw, where /j,\y denotes a scale of order mjy or m t . The Wilson coeffi- 
cients Ci(nw) pick up only small QCD corrections, which can be calculated 
within fixed-order perturbation theory. In the LL (NNLL) program, the 
matching has to be worked out at the 0(ot®) [0(cx%)] level. 

2. The evolution of these Wilson coefficients from \i = fi]y down to \i = u. b 
must then be performed with the help of the renormalization group, where 
//& is of the order of m&. As the matrix elements of the operators evaluated 
at the low scale /j& are free of large logarithms, the latter are contained in 
resummed form in the Wilson coefficients. For the LL (NNLL) calculation, 
this RGE step has to be performed using the anomalous-dimension matrix 
up to order oc l s (af). 

3. To LL (NNLL) precision, the corrections to the matrix elements of the 
operators (sj\Oi(n)\b) at the scale fj, = u.^ must be calculated to order a° 
(aj) precision. The calculation also includes bremsstrahlung corrections. 

B — > Xg'-f The error of the LL prediction of the B — > X s ~f branching frac- 
tion (J4"H I4"2" 1 I4*3" 1H3)) is dominated by a large renormalization scale dependence at 
the ±25% level, which indicates the importance of the NLL series. By conven- 
tion, the dependence on the renormalization scale fi^ is obtained by the variation 
7715/2 < //.{, < 2777,5. The three calculational steps of the NLL enterprise — Step 
1 fllSIMD, Step 2 gZHMD, and Step 3 (lJ9ll50ll5T|l52ll53j)— have been performed by 
many different groups and have been independently checked. The resulting NLL 
prediction had a small dependence on the scale /j.5 as well as on the matching 
scale /Jo below 5%. But as first observed in Ref. (54), there was a large charm 
mass-scheme dependence because the charm loop vanishes at the LL level and the 
significant charm dependence begins only at the NLL level. By varying m c /rrib 
in the conservative range 0.18 < m c /mb < 0.31, which covers both the pole 
mass value (with its numerical error) and the running mass value m c (/j c ) with 
fj, c £ [m c , mb], one finds an uncertainty of almost 10% (|55p56p . This uncertainty is 
the dominant error in the NLL prediction. The renormalization scheme for m c is 
an NNLL issue, and a complete NNLL calculation reduces this large uncertainty 
by at least a factor of two (57). This finding motivated the NNLL calculation of 
the B —7- X s 7 branching fraction. 

Following a global effort, such an NNLL calculation was recently performed 
and led to the first NNLL prediction of the B — > X s j branching fraction (|58p . 
This result is based on various highly-nontrivial perturbative calculations (159} 
[60l[6ll|62l[Mll65l[Mll66l[671[68ll69l): Within Step 1 the matching of the effective 
couplings Ci at the high-energy scale /io ~ M\y requires a three-loop calculation 
for the cases i = 7, 8 (}59~|) and a two- loop calculation for the other cases ([60]) . 
Within Step 2 the self-mixing of the four-quark operators (i = 1,...,6) and the 
self-mixing of the dipole operators [i = 7, 8) have been calculated by a three- 
loop calculation of anomalous dimensions ([61,62), and the mixing of the four- 
quark operators into the dipole operators by a four- loop calculation (j63|) . These 
two steps have established the effective couplings at the low scale fib ~ mb to 
NNLL precision. Thus, large logarithms of the form a™ +p (777.5) ^°g n ( m b/ m w), 
(p = 0, 1, 2), are resummed. Within Step 3, the calculation of the matrix elements 
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of the operators to NNLL precision, only the dominating contributions have been 
calculated or estimated by now. The dominating two- loop matrix element of 
the photonic dipole operator O7 including the bremsstrahlung contributions has 
been calculated in Refs. ([65 . 64 , 66l l67l) . The other important piece is the three- 
loop matrix elements of the four-quark operators, which has first been calculated 
within the so-called large- (3q approximation (f68|) . A calculation that goes beyond 
this approximation by employing an interpolation in the charm quark mass m c 
from m c > ntf, to the physical m c value has been presented in Ref. (|69p. In this 
interpolation the o? s j3q result (68) is assumed to be a good approximation for 
the complete a 2 result for vanishing charm mass. It is this part of the NNLL 
calculation which is still open for improvement. Indeed a complete calculation 
of the three-loop matrix elements of the four-quark operators 0\^i for vanishing 
charm mass is work in progress (|70|) and will cross-check this assumption and the 
corresponding error estimate due to the interpolation. 

Some perturbative NNLL corrections have not yet been included in the present 
NNLL estimate, but they are expected to be smaller than the current perturbative 
uncertainty of 3%: the virtual and bremsstrahlung contributions to the (07,08) 
and (08; Os) interferences at order a?, the NNLL bremsstrahlung contributions 
in the large-/3o-approximation beyond the (07,07) interference term (which are 
already available (|7ip ). the four- loop mixing of the four-quark operators into the 
operator Os (63), and the exact mass dependence of various matrix elements 
beyond the large /?o approximation (|72 y 70 tlT3"[) . 

In the present NNLL prediction (|58D , the reduction of the renormalization-scale 
dependence at the NNLL is shown in Fig. [2l The most important effect occurs 
for the charm mass MS renormalization scale ji c , which has been the main source 
of uncertainty at the NLL. The current uncertainty of ±3% due to higher-order 
[O(a^)] effects can be estimated via the NNLL curves in Fig. The reduction 
factor of the perturbative error is greater than a factor of three. The central 
value of the NNLL prediction is based on the choices //& = 2.5 GeV and \x c = 1.5 
GeV. 

At NNLL QCD accuracy subdominant perturbative electroweak two-loop cor- 
rections are also relevant and have been calculated to be —3.6% (J74"ll75"ll76[l77[) . 
They are included in the present NNLL prediction. 

B — > X s £+£~ Compared with the B — > X s ~f decay, the inclusive B — > X s £ + £~ 
decay presents a complementary and more complex test of the SM, given that 
different perturbative electroweak contributions add to the decay rate. This in- 
clusive mode is also dominated by perturbative contributions, if one eliminates 
cc resonances with the help of kinematic cuts. In the so-called perturbative q 2 - 
windows below and above the resonances, namely in the low-dilepton-mass region 
1 GeV < q 2 = rrVfe < 6 GeV as well as in the high-dilepton-mass region where 
q 2 > 14.4 GeV , theoretical predictions for the invariant mass spectrum are dom- 
inated by the perturbative contributions. A theoretical precision of order 10% is 
possible. 

Compared with the decay B — > X^, the effective Hamiltonian (Eq. [5]) con- 
tains two additional operators of 0(a era ), the semileptonic operators O9 and 
Oio- Moreover, the first large logarithm of the form log(mb/mw) already arises 
without gluons, because the operator O2 mixes into Og at oneJ.oop. It is then 
convenient to redefine the dipole and semileptonic operators via Oj = Wtt 2 /g 2 Oi, 
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Ci = g 2 /(Ai:) 2 Ci for i = 7, ..., 10. With this redefinition, one can follow the three 
calculational steps discussed above. In particular, after the reshufnings the one- 
loop mixing of the operator O2 with Og appears formally at order a s . To LL 
precision, there is only Og with a non-vanishing tree-level matrix element and a 
non- vanishing coefficient. 

It is well-known that this naive a s expansion is problematic, since the formally- 
leading 0{l/a s ) term in Cg is accidentally small and much closer in size to an 
0(1) term. Thus, also specific higher order terms in the general expansion are 
numerically important |f| 

The complete NLL contributions have been presented (j-44|.IT8[) . For the NNLL 
calculation, many components were taken over from the NLL calculation of the 
B —7- X s 7 mode. The additional components for the NNLL QCD precision have 
been calculated ([^H8l[6Tt[7^[80l[83|[8^[8^[8^[8Tt[8^[88l[87j) : Some new pieces for 
the matching to NNLL precision (Step 1) have been calculated in Ref. (|60p . To 
NNLL precision the large matching scale uncertainty of 16% at the NLL level is 
eliminated. In Step 2, the mixing of the four-quark operators into the semilep- 
tonic operator Og has been calculated (I48 p6ip . In Step 3, the four-quark matrix 
elements including the corresponding bremsstrahlung contributions have been 
calculated for the low-q 2 region in Refs. (79,80,81), bremsstrahlung contribution 
for the forward-backward asymmetry in B — > X s d. + l~ in Refs. (85,83,84), and 
the four-quark matrix elements in the high-g 2 region in Refs. (182U81ll86|) . The 
two-loop matrix element of the operator Og has been estimated using the cor- 
responding result in the decay mode B — > X u tv and also Pade approximation 
methods (88); this estimate has been further improved in Ref. (|87h . 

More recently electromagnetic corrections were calculated: NLL quantum elec- 
trodynamics (QED) two-loop corrections to the Wilson coefficients are of 0(2%) ([88 
Also, in the QED one-loop corrections to matrix elements, large collinear loga- 
rithms of the form log (m 2 /m 2 ) survive integration if only a restricted part of the 
dilepton mass spectrum is considered. These collinear logarithms add another 
contribution of order +2% in the low-q 2 region for B — > X s li + [i~ (i89i) . For the 
high-g 2 region, one finds —8% (f90|) . 

B — > Xd'Y and B — > Xd£^~£~ The perturbative QCD corrections in the 
inclusive decays B -)■ X d >y flHUEBESD and B -)■ X d i+t~ (|92l[93]) can be treated 
completely analogously to those in the corresponding b — > s modes. The effective 
Hamiltonian is the same in these processes, up to the obvious replacement of the 
s quark field by the d quark field. However, because \ u = V u bV* d for b — > d^ is 
not small with respect to At = VtbV t * d and A c = V c bV^, one must also account for 
the operators proportional to X u , namely 0" 2 in Eq. The matching conditions 
Ci{mw) and the solutions of the RGEs, which yield Oj(/i&), coincide with those 
needed for the corresponding b — > s processes (|9iy 92h. 



The B — > X s £ + £~ decay amplitude has the following structure (k — a cln /a s ): 

A — k [All + a, Anll + ce 2 s Annll + C*(a s )] with All ~ a s Anll (7) 

A strict NNLL calculation of the squared amplitude A 2 should only include terms up to order 
K 2 a 2 . However, in the numerical calculation, one also includes the term AnllAnnll of order 
k a s which are numerically important. These terms beyond the formal NNLL level are pro- 
portional to |Cr | 2 and |Cs | 2 and are scheme-independent. One can even argue that one picks 
up the dominant NNNLL QCD corrections because the missing NNNLL piece in the squared 
amplitude, namely AllAnnnll, can safely be neglected f571) . 
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2.3 Hadronic Power Corrections to Inclusive Modes 

The inclusive modes B — > X s ^f and B — > X s £ + £~ are dominated by the partonic 
contributions. Indeed, if only the leading operator in the effective Hamiltonian 
(07 for B ->■ X s 7, 9 for B ->• X s £ + is considered, the HME makes it 
possible to calculate the inclusive decay rates of a hadron containing a heavy 
quark, especially a b quark (JMIESEEEIIEEIESI)- The optical theorem relates 
the inclusive decay rate of a hadron Hb to the imaginary part of the forward 
scattering amplitude 

T(H b -> X) = —^9 (flj, | T | fT 6 ) , (8) 

2m Hb 

where the transition operator T is given by T = i J d 4 xT[H e fi(x)'H c fi(0)]. The 
insertion of a complete set of states, |X)(X|, leads to the standard formula for 
the decay rate: 

T(H b -+X) = — !— ^(2vr) 4 5 4 ( K - p f ) | (X \ U cS \ H b ) \ 2 . (9) 

lm H b x 

It is then possible to construct an OPE of the operator T, which is expressed as 
a series of local operators that are suppressed by powers of the b quark mass and 
written in terms of the b quark field: 

T[H cS n cS ] °= E —(0o + —0i + \o 2 + ...) . (10) 

m b m b ml 

This construction is based on the parton-hadron duality. The sum is performed 
over all exclusive final states and that the energy release in the decay is large with 
respect to the QCD scale, A <C m b . With the help of the heavy quark effective 
theory (HQET) (94,95), namely the new heavy quark spin-flavor symmetries that 
arise in the heavy quark limit m b — > oo, the hadronic matrix elements within the 
OPE, (H b | 0j | Hf,), can be further simplified. In this well-defined expansion, the 
free quark model is the first term in the constructed expansion in powers of l/m b 
and, therefore, is the dominant contribution. In the applications to inclusive rare 
B decays, one finds no correction of order A/m b to the free quark model approx- 
imation. The corrections to the partonic decay rate begin with l/m b only, which 
implies the rather small numerical impact of the nonperturbative corrections on 
the decay rate of inclusive modes. However, there are more subtleties to consider 
if other than the leading operators are taken into account (see below). 

B — > Xg'y These techniques can be used directly in the decay B — > X s ^ to 
single out nonperturbative corrections to the branching fraction: If one neglects 
perturbative QCD corrections and assumes that the decay B — > X s j is due to 
the leading electromagnetic dipole operator Oj alone, then the photon would 
always be emitted directly from the hard process of the b quark decay. One 
has to consider the time-ordered product TO^{x) 07(0). Using the OPE for 
TO^f{x) 07(0) and HQET methods, as discussed above, the decay width T(B — >• 
Xg'j) reads (|30U3ip (modulo higher terms in the l/m& expansion): 

r (0 7 ,0 7 ) _ aemG 2 F m 5 b 2 2 / 5^ \ NP _ 1 9 

i s->x s7 - — ^i — \VaVt8\ W( m fc) I 1 + ~^TJ> drad ~ 2 2 ' ^ ' 
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where Ai and A2 are the HQET parameters for the kinetic energy and the chro- 
me-magnetic energy, respectively. If the B — > Xg'y decay width is normalized by 
the (charmless) semileptonic one, the nonperturbative corrections of order 1/m^ 
cancel out within the ratio B{B — } X s j)/B(B — > X u lv). 

However, as noted in Ref. (I96p . there is no OPE for the inclusive decay B — > 
X s j if one considers operators beyond the leading electromagnetic dipole op- 
erator O7. Voloshin (|97l) has identified a contribution to the total decay rate 
in the interference of the 6—7-57 amplitude due to the electromagnetic dipole 
operator O-j and the charming penguin amplitude due to the current-current op- 
erator 02- This resolved photon contribution contains subprocesses in which the 
photon couples to light partons instead of connecting directly to the effective 
weak- interaction vertex ([98]) . If one treats the charm quark as heavy, then it is 
possible to expand the contribution in local operators. The first term in this ex- 
pansion may be the dominating one (96,99, 100). This nonperturbative correction 
is suppressed by X^/ml and is estimated to be of order 3% compared with the 
leading-order (perturbative) contribution to the decay rate T^^^ which arises 
from the electromagnetic operator O7: 



(02,07) 

~ +0.03. (12) 



^ L B^X s -y __ 1 O2 A2 



P^ 7 9 C 7 ml 

However, if the charm mass is assumed to scale as m 2 c ~ Am&, then the charm 
penguin contribution must be described by the matrix element of a nonlocal 
operator (|96ll99lPT00l[T0T]) . 

Recently, another example of such nonlocal matrix elements within the power- 
suppressed contributions to the decay B — > X s ^f was identified (jlOip — specifically, 
in the interference of the b — > sj transition amplitude mediated by the electro- 
magnetic dipole operator C7, where the b — )• sg amplitude is mediated by the 
chromo-magnetic dipole operator Og, followed by the fragmentation of the gluon 
into an energetic photon and a soft quark- antiquark pair. A naive dimensional 
estimate of these power corrections leads to 

Ar (e> 7 ,0s) r . 

^ET~™'w (13) 

whereas an estimate using the vacuum insertion method for the nonlocal matrix 
elements indicates an effect of —3%. 

Power corrections to the high-energy part of the B — > X s ~f photon spectrum can 
be parameterized systematically in terms of subleading shape functions. For the 
interference of the O7-O7 pair, these nonlocal operators reduce to local operators, 
if one considers the total decay rate ([102,103,104), whereas other resolved photon 
contributions to the total decay rate — such as the previously analyzed O-j-0% 
interference term (llOip — cannot be described by a local OPE. A recent systematic 
analysis ([98]) of all resolved photon contributions related to other operators in 
the weak Hamiltonian establishes this breakdown of the local OPE within the 
hadronic power corrections as a generic result. Clearly, estimating such nonlocal 
matrix elements is very difficult, and an irreducible theoretical uncertainty of 
±(4 — 5)% for the total CP averaged decay rate, defined with a photon-energy 
cut of E 7 = 1.6 GeV, remains (98). This result strongly indicates that the 
theoretical efforts for the B — > X s j mode have reached the nonperturbative 
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boundaries. The complete effect of power corrections on CP asymmetries has 
not yet been estimated. 

B — y X s £^~£~ Hadronic power corrections in the decay B — > X s £ + £~ that 
scale with 1/m 2 , l/m| (CT5T UTlMlU6lllU7l[TU8]) . and l/m 2 c (fTUU]) have also been 
considered. They can be calculated quite analogously to those in the decay 
B — > X s j. However, a systematic analysis of hadronic power corrections including 
all relevant operators has yet to be performed. Thus, an additional uncertainty 
of ±5% should be added to all theoretical predictions for this mode on the basis 
of a simple dimensional estimate. 

In the high-g 2 region, one encounters the breakdown of the HME at the end 
point of the dilepton mass spectrum: Whereas the partonic contribution vanishes, 
the 1/m 2 and l/m| corrections tend towards a nonzero value. In contrast to 
the end-point region of the photon-energy spectrum in the B — > X s j decay, no 
partial all-order resummation into a shape function is possible. However, for an 
integrated high-g 2 spectrum an effective expansion is found in inverse powers of 
m eff _ mb x (1 — ^/s m i n ) rather than nib (|109yil0|) . The expansion converges less 
rapidly, depending on the lower dilepton- mass cut s m ; n = Q r min/ rn i dSI])- 

The large theoretical uncertainties could be significantly reduced by normaliz- 
ing the B — > X s £ + £~ decay rate to the semileptonic B — > X u £v decay rate with 
the same q 2 cut (108): 

K(so) = /' u dr ' B -> y* n i f\- t ^b'-.xj.) 

Js ds Js ds 

For example, the uncertainty due to the dominating 1/mjj term would be reduced 
from 19% to 9% (1891) . 



B —¥ Xd'f The nonperturbative contributions in the decay B — > Xjj can 
be treated analogously to those in the decay B — > X s j. The power corrections 
that scale as 1/m 2 (in addition to the CKM factors) are the same for the two 
modes. Also, the systematic analysis of resolved contributions given in Ref. (|98[) 
can be applied to this case. However, the long-distance contributions from the 
intermediate u quark in the penguin loops are critical. They are suppressed in 
the B — > Xg'y mode by the CKM matrix elements. In B — > X^, there is no 
CKM suppression, and one must account for the nonperturbative contributions 
that arise from the operator 0". The contribution due to the 0\~0^ interference 
scales with A/m^ (106). However, this interference contribution vanishes in the 
total CP-averaged rate of B — > X s ^ at order A/mb ([98]) . This result applies to 
the total rate of B — > X^ as well. Other interference terms, namely the double 
resolved contributions O^-Og and 0\-0\, arise first at order 1/m 2 , as they can 
also be deduced from the results presented in Ref. ([98]) . Thus, there is no power 
correction due to the operator Of in the total rate of B — > X^ at order A/ nib, 
which implies that the CP-averaged decay rate of B — > X,fi is as theoretically 
clean as the decay rate of B — > X^. 

B — > Xd£ + £~ In the case of B — > X&£+£~ long-distance contributions due 
to u quark loops can be avoided in the low-g 2 window 1 GeV 2 < q 2 < 6 GeV 2 . 
The p and uj resonances are below, and the cc (J ftp, ip') resonances are above 
this window (92). The effect of their respective tails can be taken into account 
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within the Kriiger-Sehgal (KS) approach (see Section [23]) fllll|,lll2p . In this low- 
q 2 region, one can then treat the nonperturbative power corrections analogously 
to those in the decay B — > X s £ + £~, and one can expect a similar theoretical 
accuracy in this q 2 window. 

2.4 Nonperturbative Corrections due to Kinematical Cuts 

There are additional subtleties in inclusive modes. Kinematical cuts induce ad- 
ditional sensitivities to nonperturbative physics. 

B — > X s ~f In the measurement of the inclusive mode B — > X s j one needs cuts 
in the photon-energy spectrum to suppress the background from other B decays 

(Fig. ED. 

These shape-function effects were taken into account in the experimental anal- 
ysis, and the corresponding theoretical uncertainties due to this model depen- 
dence are reflected in the extrapolation error of the experimental results (see 
Section [5. ip . The extrapolation is done from the experimental energy cut values 
down to 1.6 GeV by use of three different theoretical schemes (114.751 111511116]) 
for averaging. 

Again constraining the analysis to the leading operator O7, a cut around 
1.6 GeV might not guarantee that a theoretical description in terms of a local 
OPE is sufficient because of the sensitivity to the scale A = mj, — 2E^ (117). 
A multiscale OPE with three short-distance scales nib, V m bA, and A has been 
proposed to connect the shape function and the local OPE region. Recently, 
such additional perturbative cutoff-related effects have been calculated to NNLL 
precision by the use of SCET methods (I118]lll9jll20p . Such perturbative ef- 
fects due to the additional scale are negligible for 1.0 GeV but of order 3% at 
1.6 GeV (118). The size of these effects at 1.6 GeV is similar to the 3% higher- 
order uncertainty in the present NNLL prediction. However, the numerical con- 
sistency of the SCET analysis has recently been questioned (112ip . Far away from 
the endpoint (Eq = 1.6 GeV), the logarithmic and nonlogarithmic terms cancel; 
the same result was presented in Ref. (1122p . Within the resummation of the 
cutoff-enhanced logarithms this feature leads to an overestimate of the 0(a^) 
terms (|12ip . Further work is needed to clarify this issue. 

There is an alternative approach to the cut effects in the photon-energy spec- 
trum that is based on dressed gluon exponentiation and on the incorporation of 
Sudakov and renormalon resummations (123. 122]) . The greater predictive power 
of this approach is related in part to the assumption that nonperturbative power 
corrections associated with the shape function follow the pattern of ambiguities 
present in the perturbative calculation (|124p . In the future, these additional per- 
turbative cut effects could be analyzed and combined together with those already 
included in the experimental average. 

B — > X s £~^~£~ In the inclusive decay B — > X s £ + £~ , the hadronic and dilepton 
invariant masses are independent kinematical quantities. A hadronic invariant- 
mass cut is imposed in the experiments (see Section [577]) . The high-dilepton- 
mass region is not affected by this cut, but in the low-dilepton mass region the 
kinematics with a jet-like X s and m\ < m&A implies the relevance of the shape 
function. A recent SCET analysis shows that by using the universality of the 
shape function, a 10 — 30% reduction in the dilepton-mass spectrum can be 



Radiative and Electroweak Penguin Decays of B Mesons 13 

accurately computed. Nevertheless, the effects of subleading shape functions 
lead to an additional uncertainty of 5% (|125yi26|) . A more recent analysis (|127p 
estimates the uncertainties due to subleading shape functions more conservatively. 
By scanning over a range of models of these functions, one finds corrections in 
the rates relative to the leading-order result to be between —10% to +10% with 
equally large uncertainties. In the future it may be possible to decrease such 
uncertainties significantly by constraining both the leading and subleading shape 
functions using the combined B — > X s ^y, B — > X U £V and B — > X s £ + £~ data (|127p . 

2.5 Charmonium Resonance Contributions 

One must also consider the on-shell cc resonances, which have to be taken out. 
Whereas in the decay B — > X s j the intermediate ip background, namely B — > 
ijjX s followed by ip — )■ X'j, is suppressed for the high-energy cut £" 7 and can be 
subtracted from the B — > X s ^ decay rate, the cc resonances show up as large 
peaks in the dilepton-invariant mass spectrum in the decay B — > X s £ + £~ . 

As discussed in Section [2T2l these resonances can be removed by making appro- 
priate kinematic cuts in the invariant mass spectrum. However, nonperturbative 
contributions away from the resonances within the perturbative windows are also 
important. In the KS approach (|llltlll2|) one absorbs factorizable long-distance 
charm rescattering effects (in which the B — > X s cc transition can be factorized 
into the product of ~sb and cc color-singlet currents) into the matrix element of 
the leading semileptonic operator 0g. Following the inclusion of nonperturbative 
corrections scaling with 1/m 2 , the KS approach avoids double-counting. For the 
integrated branching fractions one finds an increase of (1 — 2)% in the low-g 2 re- 
gion due to the KS effect, whereas in the high-g 2 region the increase is well below 
the uncertainty due to the 1/m^ corrections. As shown in Fig. O the integrated 
branching fraction is dominated by this resonance background which exceeds the 
nonresonant charm-loop contribution by two orders of magnitude. This feature 
should not be misinterpreted as a striking failure of global parton-hadron dual- 
ity ()113p . which postulates that the sum over the hadronic final states, including 
resonances, should be well approximated by a quark-level calculation ( 11281) . Cru- 
cially, the charm-resonance contributions to the decay B — > X s £ + £~ are expressed 
in terms of a phase-space integral over the absolute square of a correlator. For 
such a quantity global quark-hadron duality is not expected to hold. Neverthe- 
less, local quark-hadron duality (which, of course, also implies global duality) 
may be reestablished by resumming Coulomb-like interactions (|113p . 

2.6 Soft Collinear Effective Theory for Exclusive Decays 

The Wilson coefficients of the weak effective Hamiltonian are process-independent 
and therefore can be used directly in the description of exclusive modes. It is 
computing of the hadronic matrix elements between meson states that is diffi- 
cult in the case of exclusive modes and that limits the theoretical precision. The 
naive approach consists of writing the amplitude A ~ Cj(/Ufe)(Oj( / Ufe)) and param- 
eterizing (Oi(fj,b)} in terms of form factors. A substantial improvement can be 
obtained by using the QCDF method (32,33,34) and its field-theoretical formu- 
lation, the SCET method ([HUEgESIMllMllIID- These methods form the basis 
of the up-to-date predictions of exclusive B decays. Within this framework one 
can show that, even if the form factors were known with infinite precision, the 
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description of exclusive decays would be incomplete due to the existence of so- 
called nonfactorizable strong interaction effects that do not correspond to form 
factors. 

The QCDF and SCET methods were first systematized for nonleptonic decays 
in the heavy quark limit. In contrast to the HQET, SCET does not correspond 
to a local operator expansion. Whereas HQET is applicable to B decays if the 
energy transfer to light hadrons is small, for example to B — > D transitions at 
small recoil to the D meson, HQET is not applicable if some of the outgoing, 
light particles have momenta of order nib- If so, one faces a multi-scale problem 
that can be tackled within SCET. In this case, there are three relevant scales: (a) 
A = few x Aqcdj the soft scale set by the typical energies and momenta of the 
light degrees of freedom in the hadronic bound states; (b) m&, the hard scale set 
by both the heavy b quark mass and the energy of the final-state hadron in the 
B meson rest frame; and (c) the hard-collinear scale //he = V m bA~, which appears 
through interactions between the soft and energetic modes in the initial and 
final states. The dynamics of hard and hard-collinear modes can be described 
perturbatively in the heavy quark limit nib — > oo. Thus, SCET describes B 
decays to light hadrons with energies much larger than their masses, assuming 
that their constituents have momenta collinear to the hadron momentum. 

B — > K**y and B — >• p~y The application of the QCDF formalism to radiative 
and semileptonic decays was first proposed in Ref. (|129p . For B — > K*y, or more 
generally for B — > V7, where V is a light vector meson, the QCDF formula for 
the hadronic matrix element of each operator of the effective Hamiltonian in the 
heavy quark limit and to all orders in a s reads 

{V 1 \O i \B) = TfF B ^+ —M") dutfrMTftfau). (15) 

JO W JO 

This formula separates the process-independent nonperturbative quantities F B ~* V± 
a form factor evaluated at maximum recoil (q 2 = 0), and <ps and 4>v ± , the 
light-cone distribution amplitudes (LCDAs) for the heavy and light mesons, re- 
spectively, from the perturbatively calculable quantities T and T . The latter 
correspond to vertex and spectator corrections, respectively, and have been cal- 
culated to 0{a\) (|130U13i|ll32t[l33l) . More recently, some a 2 , terms were also 
presented (|134h . 

Light-cone wave functions of pseudo-scalar and vector mesons that enter the 
factorization formula have been studied in detail through the use of light-cone 
QCD sum rules ()135|.ll36lll37|.ll38T) . However, not much is known about the B 
meson LCDA, whose first negative moment enters the factorized amplitude at 
0(a s ). Because this moment also enters the factorized expression for the B — > 7 
form factor, it might be possible to extract its value from measurements of decays 
such as B — > yev, if the power corrections are under control. 

The QCDF formula also includes an important simplification in the form factor 
description. The B — > V form factors at large recoil have been analyzed in 
SCET (|139lll40|ll4ip and are independent of the Dirac structure of the current 
in the heavy quark limit; as a consequence, all B — > V± form factors reduce to 
a single form factor up to factorizable corrections in the heavy quark and large 
energy limits (|1421I129|) . 

Field-theoretical methods such as SCET make it possible to reach a deeper 
understanding of the QCDF approach. The various momentum regions are repre- 
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sented by different fields in the effective field theory. The hard-scattering kernels 
T 1 and T 11 can be shown to be Wilson coefficients of effective field operators. 
Using SCET one can prove the factorization formula to all orders in a s and to 
leading order in A/mb (1143[> . QCD is matched on SCET in a two-step procedure 
that separates the hard scale p ~ m b and then the hard-collinear scale p ~ \/Amj 
from the hadronic scale A. The vertex correction term T 1 involves the hard scales, 
whereas the spectator scattering term T n involves both the hard and the hard- 
collinear scales. This is why large logarithms have to be resummed (|143p . which 
can be done most efficiently in SCET. 

In principle, the field-theoretical framework of SCET allows one to go beyond 
the leading-order result in A/ nib (144). However, a breakdown of factorization is 
expected at that order (f34"|) . For example, in the analysis of B — > K*^ decays at 
subleading order, an infrared divergence is encountered in the matrix element of 
08 (|145p . In general, power corrections involve convolutions, which turn out to 
be divergent. Currently, no solution to this well-analyzed problem of end-point 
divergences within power corrections is available (|1461ll40jll47p . Thus, within 
the QCDF/SCET approach, a general, quantitative method to estimate the im- 
portant A/rrib corrections to the heavy quark limit is missing, which significantly 
limits the precision in phenomenological applications. 

Nevertheless, some very specific power corrections are still computable and are 
often numerically important. Indeed, this is the case for the annihilation and weak 
exchange amplitudes in B — > p^. The annihilation contributions also represent 
the leading contribution to isospin asymmetries (11450 . All these corrections are 
included in recent analyses of radiative exclusive decays (448,149,150). Moreover, 
the method of light-cone QCD sum rules can help provide estimates of such 
unknown subleading terms. For example, power corrections for the indirect CP 
asymmetries in B — > V7 decays have been analyzed in this manner (151). 

B — y K(*)£~\~£~ There is also a factorization formula for the exclusive semilep- 
tonic B decays, such as B — > K*£ + l~ , that are analogous to the one for the radia- 
tive decay B — > K*^/ (1301 1150]) . The simplification due to form factor relations 
is even more drastic. The hadronic form factors can be expanded in the small 
ratios A/mb and A/E, where E is the energy of the light meson. If we neglect 
corrections of order 1/mb and a s , the seven a priori independent B — > K* form 
factors reduce to two universal form factors £j_ and £11 (|142yi29j) . This reduction 
makes it possible to design interesting ratios of observables in which any soft form 
factor dependence cancels out for all dilepton masses q 2 at leading order in a s 
and A/m b (fi"52|) . 

The theoretical simplifications of the QCDF/SCET approach are restricted to 
the kinematic region in which the energy of the K* is of the order of the heavy 
quark mass; that is, q 2 <C m 2 B . Moreover, the influences of very light resonances 
below 1 GeV question the QCDF results in this region. In addition, the longitu- 
dinal amplitude in the QCDF/SCET approach generates a logarithmic divergence 
in the limit q 2 — > 0, which indicates problems in the theoretical description below 
1 GeV ()130p . Thus, the factorization formula applies well in the dilepton mass 
range 1 GeV 2 < q 2 < 6 GeV 2 . 

Clearly, the QCDF and SCET methods are also applicable to the phenomeno- 
logically important semileptonic decays such as B — > K£ + £~ (J130|.I153|) . B — > 
p£ + l~ (I150p . and B s — > (f>£ + £~ . The decay mode into a pseudoscalar is analogous 
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to the decay mode into a longitudinal vector meson. 

3 EXPERIMENTAL TECHNIQUES 

The T(45) resonance produced by the e + e~ collision at the B factories provides 
a clean sample of B°B° and B + B~ meson pairs as well as strong kinematical 
constraints that are otherwise unavailable, particularly at hadron colliders. The 
main background is from continuum light quark pair production (e + e~ — >• qq, 
q = u,d, s, c), which has a cross section only three times larger than that of BB 
production. Radiative and electroweak penguin B decays are efficiently measured 
at the B factories thanks to their clear signatures: a high-energy photon and a 
lepton pair, respectively. 

3.1 Exclusive B Reconstruction 

A B meson decaying into an exclusive final state is reconstructed by measuring all 
long-lived decay products {n^, K^, e^, p^ and 7), selecting intermediate states 
of certain invariant masses, and calculating two standard variables: the beam- 
energy constrained mass Mb c = \/ s /4 ~ \Pb\ 2 ( a ls° referred to as the beam-energy 
substituted mass, Mes) and the energy difference AE = E* B — y/s/2. Here, \/ r s/2 
is the beam energy, and p B and E B are the momentum and energy, respectively, 
of the reconstructed B meson candidate in the T(4S) rest frame. 

Mb c has a peak at the B meson mass and AE has a peak at zero (Fig. H]). 
The resolution of Mb c is significantly better than that of AE, as the former is 
dominated by the spread of the beam energy, whereas the latter is dominated 
by the detector resolution. The AE variable is sensitive to misreconstructed 
background events, whereas M^ c has little separation power for them. When a 
kaon is misidentified as a pion, AE shifts by approximately 50 MeV, and when a 
low momentum pion is added or missed, AE shifts by more than the pion mass. 
An exclusive b — > d~f final state is thus separated from a similar b — > 57 state 
with AE, but this separation is marginal due to the photon energy resolution of 
the electromagnetic calorimeter. Therefore, pion to kaon separation is crucial for 
the measurement of the suppressed b — > dj processes. 

Background events due to random combinations of particles are also reduced 
by correctly identifying particle species. For the b — > s£ + £~ processes, electrons 
and muons are almost completely separated from the more abundant hadrons. 
In addition, various techniques based on the event topology can be applied to 
suppress the background from continuum qq events. 

3.2 Inclusive Measurement Techniques 

A fully inclusive measurement of B — > X s ^, in which the system recoiling against 
the emitted photon is not reconstructed, has been performed at the B factories 
thanks to the clean environment. The dominant background photon sources are 
(a) the copiously produced ir° — > 77 decays, (b) 77 — > 77 to a lesser extent, and 
(c) other secondary and initial-state radiation photons in continuum qq events. 
These contributions can be safely subtracted because they are measured in events 
taken 60 or 40 MeV below the T(4S) resonance (i.e., off resonance). Here, small 
corrections due to the center-of-mass energy difference are applied to the pro- 
duction cross section and the reconstruction efficiency. To avoid sacrificing the 
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B decay sample for other studies, the size of the off-resonance data sample is 
only ~10% of the on-resonance sample from both Belle and BaBar and is the 
dominant source of statistical and systematic errors (CLEO collected one-third 
of the sample as off-resonance) . The second severe background source arises from 
similar secondary photons from B decays. These contributions are subtracted 
from the expected photon spectrum on the basis of measured 7r° and 7] spectra 
from B decays, various control samples, or Monte Carlo simulation. 

An alternative technique is to measure as many exclusive modes as possible and 
then calculate their sum (i.e., the sum-of-exclusive method). Exclusive branching 
fractions measured to date do not saturate the inclusive process, but one can still 
infer the total branching fraction by estimating the fraction of unmeasured modes 
of typically ~45 % (or ~30% if K L modes are accounted for by corresponding Kg 
modes) using simulated hadronization processes. In the simulation, a light quark 
pair is generated according to the SM mass spectrum and final-state hadrons are 
produced by the PYTHIA program (|155p . This method also provides direct in- 
formation about the B meson. For example, the B meson momentum defines the 
B meson rest frame, and charge and flavor information allows CP- and isospin- 
asymmetry measurements. So far, the sum-of-exclusive method is the only way 
to perform inclusive measurements of B — > X s l + £~ and B — > X^j decays. 

Another potentially definitive method is the so-called B-reco technique, in 
which the other B meson is fully reconstructed, thereby allowing the target B 
decay to be measured in a very clean environment. The efficiency is as low as a 
fraction of a percent, and will be more important in future experiments. 

4 PRESENT THEORETICAL PREDICTIONS 

Theoretical predictions have significantly improved in the past decade along with 
the development of the theoretical tools. There have also been improvements in 
the relevant experimental input quantities, as discussed below. 



4.1 Inclusive Penguin Decays 

The inclusive radiative and electroweak penguin modes offer theoretically clean 
observables because nonperturbative corrections are small and well under control. 
This assessment also applies to the branching fraction of B — > X^ mode as 
discussed in Section I 



Inclusive B — > X s j The stringent bounds obtained from B — > X s j on 
various nonstandard scenarios are a clear example of the importance of clean 
FCNC observables for discriminating NP models. 

The branching fraction for B — > X q ^f (q = s, d) can be parameterized as 



B(B -► X ql ) Ei>Eo = B{B -► X c ev) 
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P(E ) + N(E ) , (16) 



where a em = cQ she11 (74), C = \V ub \ 2 /\V cb \ 2 x T[B -* X c ev)/T\B -► X u ev) and 
P{Eq) and N(Eq) denote the perturbative and nonperturbative contributions, 
respectively. The latter are normalized to the charmless semileptonic rate to 
separate the charm dependence. 
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The first NNLL prediction, which is based on the perturbative calculations 
discussed in Section 12,21 and on the analyses of nonperturbative corrections pre- 
sented in Sections l2.3l and [2T¥| for a photon-energy cut E^ > 1.6 GeV ([58]) . reads 
as: 

B(B -»■ X s7 )nnll = (3.15 ± 0.23) x 1(T 4 . (17) 

The overall uncertainty consists of nonperturbative (5%), parametric (3%), per- 
turbative (scale) (3%) and m c -interpolation ambiguity (3%), which are added 
in quadrature. An additional scheme dependence in the determination of the 
prefactor C has been found ([156); it is within the perturbative uncertainty of 
3% (fT2T|) . 

Thus, the SM prediction and the experimental average (see Section 15. ip are 
consistent at the 1.2cr level. This finding implies very stringent constraints on 
NP models, such as (a) the bound on the charged Higgs mass in the two-Higgs 
doublet model (JI52II58D ( M H+ > 295GeV at 95% CL) (J58D and (b) the bound 
on the inverse compactification radius of the minimal universal extra dimension 
model (l/R > 600GeV at 95% CL) (p.59). In both cases, the bounds are much 
stronger than those derived from other measurements. Constraints within various 
supersymmetric extensions have been analyzed in Refs. (|160U161 ,162.163.164. T65| 
I166pi67pi68p (for overviews see (I16 pi69p ). Bounds on the little Higgs model with 
T-parity have also been presented (|170p . Finally, model- independent analyses 
in the effective field theory approach without (|17ip and with the assumption of 
minimal flavor violation (|172U173]l also show the strong constraining power of the 
B —7- X s 7 branching fraction. 

Inclusive B — > Xd"f The theoretical predictions for the branching fraction 
B{B — > Xd'y) for photon energies E-y > 1.6 GeV read as (55,56): 

B(B -► X dl ) = (l.38 ±{jif ^ ± 0.15 C KM ± 0.09 paran , ± 0.05 SC£ j c ) x 10~ 5 , (18) 



m b 



and 

BjB -j. X dl ) _ Nv ., +()l , 
B(B -> X sl ) 



(3.82 i^j^ ±0.42 C km ±0.08 param . ±0.15 scalc ) x lO" 2 . (19) 



These predictions are of NLL order. They are fully consistent with previous 
results (|9ip . A good part of the uncertainties cancel out in the ratio. The errors 
are dominated by CKM uncertainties, and thus the measurement of B{B — > X^) 
constrains the CKM parameters. This measurement is also of specific interest 
with respect to NP, because its CKM suppression by the factor |14d| 2 /|Vt s | 2 in 
the SM may not hold in extended models. 

Direct CP Asymmetry Other important observables are the direct CP 
asymmetries (q = s,d), whose sign is always defined in terms of b — b, or 

A (71 : x ,_ nB^X ql )-T(B^X- ql ) 

As first noted in Ref. (|174p . the SM predictions are almost independent from the 
photon energy cut-off and, for E^ > 1.6 GeV, read as (|55|. I56 |) 

Acp(B -> X sl ) = (0.44 t° f mc ± 0.03ckm ±Ul , ) x 10" 2 , (21) 



scale 
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and 

Acp(B -> X dl ) = (-10.2 tl A 7 mc ± I.Ockm tli , ) x 10- 2 . (22) 

V — £ scale / 

m b 

The two CP asymmetries are connected by the relative CKM factor A 2 [(1 — p) 2 + 
rj 2 ]. The small SM prediction for the CP asymmetry in the decay B — > X s ~/ is a 
result of three suppression factors: (a) a s to have a strong phase; (b) CKM 
suppression of order A 2 ; and (c) GIM suppression of order (m c /mb) 2 , which 
reflects that in the limit m c = m u , any CP asymmetry in the SM would vanish. 
On the basis of CKM unitarity, one can derive the following [7-spin relation 
between the un- normalized CP asymmetries (1175h : 

T(B -»■ X sl ) -T(B -»■ X- sl )] + [V(B -»■ X dl ) -T{B -»■ X 2l )] = (23) 

[/-spin breaking effects can be estimated within the HME (even beyond the par- 
tonic level), so one arrives at the following prediction for the total (or untagged) 
B ->■ X s+d -f asymmetry (1176111770 : 

\AB(B ->■ X sl ) + AB(B -> X d ~/)\ ~ 1 • 10" 9 . (24) 

Because this null test is based on the CKM unitarity, it represents a clear test 
for new CP phases beyond the CKM phase (|1761I177|) . NP sensitivities of direct 
CP asymmetries have been analyzed (J174|.l55j) . 

Inclusive B — > X s £~^~£~ The decay B — y X s £ + £~ is particularly attractive 
because it offers several kinematic observables. The angular decomposition of the 
decay rate provides three independent observables, Ht, Ha and Hl, from which 
one can extract the short-distance electroweak Wilson coefficients that test for 
NP (fTTHD : 

^- = -[(1 + z 2 )H T (q 2 ) + 2(1 - z 2 )H L (q 2 ) + 2zH A (q 2 )} . (25) 

Here z = cos^, 0£ is the angle between the negatively charged lepton and the 
B meson in the center-of-mass frame of the dilepton system, and q 2 is the dilep- 
ton mass squared. Ha is equivalent to the forward-backward asymmetry, and 
the dilepton-mass spectrum is given by Ht + Hl- The observables depend on 
the Wilson coefficients CV, Cg and C\q in the SM. The present measurements 
of the B —)■ X s £ + £~ already favor the SM-sign of the coefficient CV, which is 
undetermined by the B — > X s ^ mode (|179p . 

As discussed above, these observables are dominated by perturbative contribu- 
tions in the perturbative low- and high-g 2 windows which are below (1 GeV < 
q 2 < 6 GeV 2 ), and above (q 2 > 14.4 GeV 2 ) the cc resonances, respectively. The 
present predictions are based on the perturbative calculations to NNLL preci- 
sion in QCD and to NLL precision in QED (see Section [2.20 . For the branching 
fraction in the low-q 2 region one arrives at (j89[) 

K(B^X£+f-\ _ f (1-59 ±0.11) xlO- 6 (l = n) , . 

B(B -+ X s £ £ ) low - | (1 ^ ± Q n) x 1Q _ 6 (£ = eh (26) 

and for the high-q 2 region, one arrives at (|90p 



E(E^X£+£-\ _ / 2-40 x 10- 7 x (1^) {£ = p) 

B(B -> X s £ £ ) high - | 2 og x iQ _ 7 x (il o.32 } {l = e y (27) 
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As suggested in Ref. ([108), normalizing the B — > X s £ + £~ decay rate in the high- 
er 2 region to the semileptonic B — > X u £v decay rate with the same q 2 cut (Eq. [14"]) , 
significantly reduces the nonperturbative uncertainties ([90]) : 

^o-14.4GeV)-| 194xl0 _3 x(1±016) {£ = eh (28) 

The value of q^ for which the forward-backward asymmetry vanishes, 

(a 2 )\X £+£-]- I (3-50 ± 0.12) GeV 2 (£ = M ) , . 

(q )[X s £ I \~\ {338 ± QU) QeV 2 {£ = eh (29) 

is one of the most precise predictions in flavor physics and also determines the 
relative sign and magnitude of the coefficients CV and Cg (|90p . However, unknown 
subleading nonperturbative corrections of order 0(a s A/n r n J ), which are estimated 
to give an additional uncertainty of order 5%, have to be added in all observables 
of the B ->■ X s t + f mode (see Section E3|) . 

In all predictions, it is assumed that there is no cut in the hadronic mass 
region (see Section I2.4p . Furthermore, after including the NLL QED matrix 
elements, the electron and muon channels receive different contributions due to 
terms involving ln(m 2 /mf) (see Section 12. 2[) . This is the only source of the 
difference between these two channels. All collinear photons are assumed to be 
included in the X s system, and the dilepton invariant mass does not contain any 
photons; in other words, q 2 = (pi+ +Pe~) 2 - Present experimental settings at 
the B factories are different, and therefore the theoretical predictions have to be 
modified fl87|). 

This difference in the settings also means that deviations from the SM predic- 
tion (R^ = 1) in the muon-electron ratio 

R Xs = T(B -»■ X s n+u.-) [qlql] /T{B -»■ X s e + e-) [qUl] (30) 

can result from a different treatment of collinear photons in the two modes. This 
ratio is interesting because it is sensitive to the neutral Higgs boson of two-Higgs- 
doublet models at large tan (3 (]180lll81j) , which is also valid in corresponding ratios 
R K (*) of exclusive modes: In the SM, one finds Rk = 1, as well as Rk* = 0.75 
when integrated over all q 2 , including M e + e - < 2m^. 

4.2 Exclusive Penguin Decays 

The exclusive penguin modes offer a larger variety of experimentally accessible 
observables than do the inclusive ones, but the nonperturbative uncertainties in 
the theoretical predictions are in general sizable. 

B — > K*~f and B —t p*y The large hadronic uncertainties, which arise from 
the nonperturbative input of the QCDF formula and from our limited knowledge 
of power corrections, do not allow precise predictions of the branching fractions of 
exclusive modes. However, within ratios of exclusive modes such as asymmetries, 
parts of the uncertainties cancel out and one may hope for higher precision. 

The ratio i?th(P7 / 'K* 7) [and similarly R t h(^7 / K* 7)] is given by I|148|I149U150| 
[151]). 



D ( IV* \ £th(£ ~± Pi) Q 



Vm 



V t 



Is 



2 (M 2 -m 2 ) 3 „ 

( M B- m K*r 

(31) 



and 



Radiative and Electroweak Penguin Decays of B Mesons 21 

where m p is the mass of the p meson; £ is the ratio of the transition form factors, 
£ = T^(0)/Tf-*(0); and S p = 1 and 1/2 for the p^ and p° mesons, respectively. 
The quantity (1 + AR) entails the explicit 0(a s ) corrections as well as the power- 
suppressed contributions. These functions also depend on CKM parameters, 
namely 4>2 = a = arg(-V cb V* b /VtdV* b ) and R ut = \V ud V* b /V t dV t l\, and one finds 
numerically (150) that 

AR(p ± /K* ± ) = {l - 2R ut cos ( i )2 [0.24±g;il] + R 2 ut [0.07±g;^]} , (32) 

AR(p°/K*°) = {l- 2R ut co S <p 2 [-0.06±8:8|] +^t [0.02±g;g=]} . (33) 

These results are consistent with the predictions given in Refs. (|148|,ll49|ll5ip . 
Obviously, the neutral mode is better suited for the determination of |V^/Vj s | 
than is the charged mode, in which the function AR is dominated by the weak- 
annihilation contribution, which leads to a larger error. The most recent deter- 
mination of the ratio £ = T^(0)/T{^*(0) within the light-cone QCD sum rule 
approach (|182p . 1/Q = 1.17 ± 0.09, leads to the determination of |T4d/T4 s | via 
Eq. (|3ip (see Section 15. 4p . However, the experimental data on the branching 
fractions of B — > K*^ and B — >• pj calls for a larger error on £, if one assumes 
no large power corrections beyond the known annihilation terms (|150p . 

Isospin Asymmetry in Radiative Decays Another important observable 
is the isospin breaking ratio given by 

a (r^k*\ r(ff ^A* 7 )-r(£+^A*+ 7 ) 

A 0+ (B -> K i) = r{B0 - R ^ - r{B+ - ^ + ^ , (34) 

where the partial decay rates are CP-averaged. In the SM spectator-dependent 
effects enter only at the order A/ mb, whereas isospin-breaking in the form fac- 
tors is expected to be a negligible effect. Therefore, the SM prediction is as 
small as 0(5%) ljl45lll48lH5giU5DlU5Hl . The ratio is especially sensitive to NP 
effects in the penguin sector, namely to the ratio of the two effective couplings 
Cq/C-j. The analogous isospin ratio in the p sector strongly depends on CKM 
parameters (150): 



T{B+ -» p+ 7 ) _ _ ,_. +5.4, +5.8, s . n -2 

2T(B° -)• p°j) ' l) -4.2 1 CKM -5.6 1 had i 



A(P7) = :r Dn "o! - 1 = (-4.6±|i| CKM ±|;l| liad ) x 10- 2 . (35) 



The hadronic error is due mainly to the weak-annihilation contribution, to which 
a 50% error is assigned. 

CP asymmetries in Radiative Decays In the CP asymmetries, the un- 
certainties due to form factors cancel out to a large extent. But both the scale 
dependence and the dependence on the charm quark of the next-to-leading-order 
predictions are rather large because the CP asymmetries arise at 0{a s ) only. Al- 
though the direct CP asymmetry in B — > K*^f is doubly Cabibbo suppressed and 
expected to be very small within the QCDF/SCET approach, one finds 0{— 10%) 
predictions for the direct CP asymmetries in the B -> p 7 mode (113ipi48lll50p . 
Because the weak-annihilation contribution does not contribute significantly here, 
the neutral and charged modes are of similar sizes (jl50[> : 

Acp(B° -> A) = (-10.4i^| CKM t^| had ) % (36) 
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and 

Acp(B- -> p~i) = (-10.7^; 5 | CKM ±! 6 7 | had )%. (37) 

Finally, we reiterate that all predictions of exclusive observables within the QCDF/ 
SCET approach may receive further uncertainties due to the unknown power 
corrections. This possibility might be especially important in the case of CP 
asymmetries. 

The time-dependent CP asymmetry is given by two parameters, Sep and Acp'- 

A C p{B ->■ /; At) = S C p sin(AmAi) + A C p cos(AmAt), (38) 

where Aqp represents the size of the direct CP asymmetry discussed aboveo In 
hadronic decay modes such as B — > J/ipKg, a large value of Sep due to the angle 
<pi = j3 = — &Yg(V t dV t * b /V uc iV* b ) of the unitarity triangle has been established, and 
a similarly large CP asymmetry is expected for hadronic penguin decays. This 
asymmetry is suppressed in radiative penguin decays because the photon helicities 
are opposite between those from B and B° decays under the left-handed current 
of SM weak decays, and they do not interfere in the limit of massless quarks. 
This finding implies a suppression factor of m s /rrib in the leading contribution to 
Sep that is induced by the electromagnetic dipole operator O7: 

5 SM = _ gin 2(j)i ™£ [2 + 0{ag)] + 5 SM, S79 (39) 

m b 

However, there are also additional contributions, S ,a ^ 9 induced by the process 
b — > s'jg via operators other than O7 (183 11841) . These corrections are not helicity- 
suppressed but are power-suppressed. A conservative dimensional estimate of the 
contribution from a nonlocal SCET operator series leads to |5 ,S79 | ~ 0.06 (|183^ 
184), whereas within a QCD sum rule calculation, the contribution due to soft- 
gluon emission is estimated to be 5 SM ' S79 = -0.005 ± 0.01 (JTHKIITBT]) which 
leads to Sjj$ = -0.022 ± 0.015 ig_ 01 @ The QCD sum rule estimates of power 
corrections, namely long-distance contributions that arise from photon and soft- 
gluon emission from quark loops (|15ip . lead to analogous results for the other 
radiative decay modes, such as B — > /r/ (I15ip . If a large value of Sep beyond 
the SM prediction is observed, it will signal a new right-handed current beyond 
the SM. 

B — > K*£~^~£~ The isospin asymmetry in the mode B — > K*C + l~ , as in the 
radiative mode, is a subleading A/mb effect, but the dominant isospin-breaking 
effects can be calculated perturbatively, whereas other A/m& corrections are sim- 
ply estimated. Thus, the exact uncertainty is difficult to estimate due to unknown 
power corrections, but the observable may still be useful in the NP search because 
of its high sensitivity to specific Wilson coefficients (186). 

The decay B° — > K*°£ + £~ (with K*° — > K~tt + on the mass shell) is completely 
described by four independent kinematic variables: the lepton-pair invariant mass 
squared, q 2 , and the three angles 0£, 9k, and (j) (for their precise definitions, see 



4 The symbol Cap = —Acp is also often used. 

5 This does not necessarily contradict a larger time-dependent CP asymmetry of approxi- 
mately 10% within the inclusive mode found in Ref. (|183[) . because the SCET estimate 1|183II184[) 
shows that the expansion parameter is A/Q. Here Q is the kinetic energy of the hadronic part. 
There is no contribution at leading order. Thus, the effect is expected to be larger for larger 
invariant hadronic mass. The K* mode must have the smallest effect, below the average 10%. 
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Ref. (|152pi87l) ). Summing over the spins of the final particles, the differential 
decay distribution can be written as (|188yi89|) 

~ -—I(q 2 ,e i ,e K ,<P) S me i sm6 K . (40) 



dq 2 dOi d9 K d<f> 32?r 
By integrating two of the angles, one finds 
dF' 3V 



dB K 4 



sin 9 K 2F l cos^ 6 k + (1 - F L ) sin^ 9 K , (41) 



and 



— = T' ( ~F L sin 2 9 t + -(1 - F L )(1 + cos 2 0,) + A FB cos^J sin0 € . (42) 

The observables appear linearly in the expressions so the fits can be performed 
on data binned in q 2 . The fraction of longitudinal polarization Fl from the kaon 
angular distribution and the forward-backward asymmetry Afb from the lepton 
angular distribution are accessible this way. The latter observable is defined as 
follows (Qn is defined below Eq. [25]) : 

1 / C 1 d 2 T r° d 2 T \ 

The hadronic uncertainties of these two differential observables are large. How- 
ever, the value of the dilepton invariant mass q 2 , for which the differential forward- 
backward asymmetry vanishes, can be predicted in quite a clean way. In the 
QCDF approach at leading order in A/m^, the value of q 2 , is free from hadronic 
uncertainties at order a®. A dependence on the soft form factor and on the 
light-cone wave functions of the B and K* mesons appears only at order a\. At 
next-to-leading order one finds ([150) : 

q^[K*°e + r] = 4.36±g;|? GeV 2 , q 2 [K* + e + r] = 4.15J£|? GeV 2 . (44) 

The small difference is due to isospin-breaking power corrections. However, an un- 
certainty due to unknown power corrections should be readded to the theoretical 
error bars. The zero is highly sensitive to the ratio of the two Wilson coefficients 
C7 and Cg. Thus, such a measurement would have a huge phenomenological 
impact. 

In the near future, a full angular analysis based on the four-fold differential 
decay rate in Eq. [JO] will become possible. Such rich information would allow 
for the design of observables with specific NP sensitivity and reduced hadronic 
uncertainties (|152]I187]) . These observables would be constructed in such a way 
that the soft form factor dependence would cancel out at leading order for all 
dilepton masses, and they would have much higher sensitivity to new right-handed 
currents than would observables that are already accessible via the projection 
fits (|189|ll52|ll87j) . In these optimized observables, the unknown A/m^ corrections 
would be the source of the largest uncertainty. Further detailed NP analyses of 
such angular observables have been presented in Refs. (|19iyi90|) . A full angular 
analysis provides high sensitivity to various Wilson coefficients, but the sensitivity 
to new weak phases is restricted (|192yi87|) . 
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5 PRESENT EXPERIMENTAL RESULTS 

The huge samples of B meson decays collected by Belle and BaBar have made 
it possible to fully explore the radiative penguin decays b — > 57 and b — > dy, as 
well as the electroweak penguin decays b — > s£ + £~ . 

5.1 Inclusive B — V X a q Branching Fraction 

An experimental challenge is how to lower the minimum photon energy to 1.6 GeV 
(see Section 12.4ft . Before the construction of the B factories, minimum photon 
energy of 2.0 GeV was required in the measurement by CLEO (193). BaBar 
has a minimum photon-energy requirement of 1.9 GeV based on 89 million BB 
pairs (|194p . whereas Belle first reported the result of 1.8 GeV with 152 million 
BB (|195j) . Belle recently lowered the limit to 1.7 GeV by using 657 million BB 
pairs (196) (Fig. [5]). Belle measured the branching fraction to be B(B — > X s j) = 
(345 ± 15 ± 40) x 10~ 6 for E 1 > 1.7 GeV, whereas BaBar measured it to be 
(367 ± 29 ± 34 ± 29) x 10~ 6 for £ 7 > 1.9 GeV. 

The sum-of-exclusive technique and the B-reco technique have been used by 
BaBar with a minimum photon-energy requirement of 1.9 GeV (|197tll98"j) . which 
corresponds to a maximum recoil mass requirement of 2.8 GeV. Belle also made 
a sum-of-exclusive measurement using a very early data set (199). 

To calculate the average branching fraction based on the same phase space 
and to compare it with theory predictions, the Heavy Flavor Averaging Group 
(HFAG) (200) has made an extrapolation of the branching fraction to the same 
minimum photon energy of 1.6 GeV (|116j) . The extrapolation factor is 0.985 ± 
0.004 (0.894 ± 0.016) for 1.7 (2.0) GeV. The average thereby obtained is 

B{B ^ X s -f) = (352 ±23 ±9) x 10~ 6 , (45) 

where the first error is statistical and systematic combined, and the second is due 
to the extrapolation. The result is in agreement with the SM prediction given in 
Eq. [T71 and it provides stringent constraints on NP, as discussed in Section 14.11 

5.2 Exclusive Measurements of b — > 57 Processes 

The recoil system of B — > X s j below 1.1 GeV is dominated by the K* resonance, 
as the spin-0 state is forbidden. Above 1.1 GeV, X s is a mixture of various reso- 
nant and nonresonant states and therefore can usually be modeled as a continuum 
spectrum in the inclusive B — > X s ^f analysis. The K* signal makes it possible to 
use the channel for various studies, as discussed below. The B — > K*^f branching 
fractions have been measured precisely by Belle (|20ip and BaBar (I154D and have 
been averaged by HFAG to be 

B(B° -)■ K*°j) = (43.3 ± 1.5) x lO" 6 , . 

B{B+ -> K*+j) = (42.1 ± 1.8) x 10~ 6 , ^ > 

which corresponds to approximately 12% of the total B — > X s ^ branching frac- 
tion. The SM predictions for the branching fraction have been calculated by 
many groups. They are consistent with the measured values but have very large 
errors of 30% to 50%, which arise mainly from the uncertainty of the B — > K* 
form factor p30U13ip . 
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The resonant structure of the high mass X s system has also been explored. So 
far, B -> K| (1430)7 d2Q21[203D and B -»• K x (1270)7 (gOU) have been measured, 
but other decay channels such as B — > K\ (1400)7 seem to have small branching 
fractions and have not yet been observed. In addition, many multi-body final 
states have been measured, such as B — > Kmr'y (|205P202"1) (including B — > K*tt^ 
and B -> K/rf), B -»• Kcjry J2JEHMD, B -> ^777 ([208l[2MD, 5 -»• #7/7 (J2T0|) . 
and 5 -*■ Ap-y pTTj) . 

5.3 CP and Isospin Asymmetries in 6 — > 57 Processes 

The measurement of the inclusive direct CP asymmetry (Eq. [20]) was performed 
by use of the sum-of-exclusive method to tag the flavor of the B candidate. For 
B° (B°), only the self-tagging modes with a K + (K~) were used. The measured 
asymmetry was corrected for a small dilution due to the doubly misidentified pair 
of a charged kaon and a pion. The results, based on 152 and 383 million BB 
samples by Belle and BaBar, are 0.002±0.050±0.030 (|2"T2"j) and -0.011±0.030± 
0.014 (I213p . respectively, and have been averaged by HFAG to be 

A C p(B -> X sl ) = -0.012 ± 0.028. (47) 

This is consistent with null asymmetry, and the size of the error is still much larger 
than the SM precision (Eq. 12 1|) . Assuming the systematic error can be reduced 
along with the statistical error, a data set two orders of magnitude larger would 
be more sensitive to NP, although still insufficient to measure the small Acp 
predicted by the SM. 

In the exclusive B — > K*"f channel, the direct CP asymmetry is also small but 
has less-understood theoretical uncertainties. Experimentally it can be measured 
more precisely. The current HFAG average is 

A C p{B ->• K*j) = -0.003 ± 0.017, (48) 

which is also consistent with null asymmetry. 

In the fully inclusive measurement, flavor information is not available for the 
signal side, but it can be obtained from the charge of the lepton in the event if the 
other B decays into a semileptonic final state. In this case, it is not possible to 
discriminate B — > Xd'y from B — > X s j and the measured asymmetry corresponds 
to a combined one. This combined asymmetry (note the different normalization 
in comparison with Eq. [24"|) has been measured by BaBar to be 

Acp{B -> A s+d7 ) = -0.110 ± 0.115 ± 0.017, (49) 

which is consistent with null asymmetry but has a much larger error than do the 
other two asymmetry measurements. 

The measurement of the isospin asymmetry (Eq. [34"|) is another way to utilize 
the reconstructed B — > K*j events. Here, the measured branching fractions are 
corrected by the lifetime ratio: t b +/t b o = 1.071 ± 0.009 (|214p . Usually B decay 
branching fractions are quoted based on the assumption of B(T(AS) — > B + B~) = 
B(T(4S) —> B°B°) = 0.5, but the isospin-asymmetry measurement is already 
precise enough to be affected by the difference in these branching fractions. Belle 
measures Aq + (B — > K*^) = 0.012±0.044±0.026 without this correction, whereas 
BaBar measures A 0+ (B ->■ K*-f) = 0.066 ±_0.021 ± 0.022 using B(T(4S) -»• 
B+B~) = 0.516 ± 0.006 and B(T(4S) ->■ B°B°) = 0.484 ± 0.006 (|214"|) . After 
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scaling the Belle result and including the CLEO result using the aforementioned 
lifetime and production ratios, the naive world average is 

A 0+ (B ^ K*j) = 0.062 ±0.027. (50) 

This average is in agreement with the SM expectation. 

A similar isospin asymmetry can be also measured for the inclusive B —> X s ~f 
decay by use of the sum-of-exclusive method. The result by BaBar is 

A 0+ (B ->• X s7 ) = -0.006 ± 0.058 ± 0.009 ± 0.024, (51) 

which is consistent with null asymmetry but is not yet as precise as that for 
B -> K*j. 

The measurement of the time-dependent CP asymmetry (Eq. 1381) for b — > 
s-f faces two experimental challenges. First, the modes and statistics that can 
be used for time-dependent CP asymmetry measurements are rather limited. 
Although the B° — > K* 0r y branching fraction is not very small, only 1/9 of the 
events that decay into the Kg(— > 7r + 7r _ )7r°7 final state can be used. Second, 
the B meson decay vertex position has to be extrapolated from the displaced 
Kg —7- 7r + 7r~ vertex and the Kg momentum vector. Therefore, the Kg decays 
inside the vertex detector volume (55% in Belle, 68% in BaBar) and the resulting 
vertex resolution is somewhat degraded. 

Because any B — > P°Q 0, y final states [where P° and Qp are CP eigenstates (|215h ] 
can be used, the B — > K s tt ^ events with M K o n o up to 1.8 GeV including 
-ftT|(1430) were measured by both Belle (216) and BaBar (I217p and have been 
averaged by HFAG as 

S C p{B -» Jf§7T°7) = -0.15 ± 0.20, (52) 

in which the B° -)• K* ^ contribution gives S C p(B -> K*j) = -0.16 ± 0.22. 

As additional channels, BaBar has measured Scp(B — > K s r]j) = — 0.18_ ' 46 ± 
0.12 p)9]) . and Belle has measured S C p{B -»• K s p°j) = 0.11 ± 0.33 jftgg (J23SI). 
The latter is slightly diluted by the B — > K* + ir~ , y events, by a factor which 
was measured to be 0.83l ' 3, but is free from the restriction of K$ vertexing 
and has a statistical error comparable in size to that of the B — > Kgir°j mode. 
Currently all results are compatible with null asymmetry with errors that are still 
not small enough to provide nontrivial constraints on right-handed currents, but 
this observable will be one of the best ways to search for NP in future experiments. 

5.4 Measurements of b — > dj Processes 

Three exclusive b — > d'y decay modes are considered to be the easiest modes to 
study the b — > d'y process: B + — > p + 7, B° — > /9°7 and B° — > wy. Although these 
modes have been searched for since the beginning of Belle and BaBar, only in 
the later stage of these experiments were measurements of the B — > p7 modes 
established. This is partly because of the large B — > K*^ background and partly 
due to the huge continuum background, which is more severe for modes without a 
kaon in the final state. Therefore, large statistics and good particle identification 
are essential; Belle has the advantage in the former, whereas BaBar leads in the 
latter. 

To gain statistics, these three modes have been combined by assuming their 
naive quark contents, through the use of T(B + — > p + ^i) = 2T(B° — > p°j) = 
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2T(B° -)■ urf). In the latest measurements by both Belle ([2151) and BaBar ([220]) . 
the 5° — > /3°7 mode was measured with more than 5<r significance and B + — > p + ^f 
with more than 3d significance, whereas B° — > wy remains unestablished with 
significance less than 3a. Using the symbol B — > (p, 07)7 for the combined results 
that are adjusted for the B + — > p +r y mode, the averaged branching fraction by 
HFAG becomes 

B(B -> (p,w)7) = (l-30l° i 8 9 ) x 10- 6 . (53) 

The results are consistent with the SM predictions. However, these predictions 
are affected by form factor uncertainties and do not have effective prediction 
power for NP. 

A more effective way to use these results is to combine them with the B — > K*^ 
measurements to determine iV^/Vy. Using Eq. [3TI Belle and BaBar reported 
the value of \V U /V ts \ to be 0.195 ±£019 ±0.015 and 0.233 jftra! ±omi, respectively, 
where the errors are experimental and theoretical. These results from the penguin 
diagrams are in agreement with the determination from the box diagrams using 
the ratio of the B° and B® mixing parameters Am,i/Am s , where Am^ is measured 
at the B factories and Am s at the Tevatron. The results are also in agreement 
with the more indirect determination by the unitarity triangle fit from other 
observables. This is a nontrivial test of the CKM scheme. However, although 
the experimental errors are still larger than the theoretical errors, the size of the 
theoretical error is unlikely to be reduced. 

A possible way to improve this situation utilizes the inclusive B — > X^ mea- 
surement with the sum-of-exclusive method. BaBar has reconstructed the X^n 
system in seven final states (7r + 7r~7, 7r + 7r°7, 7r + 7r~7r + 7, 7r + 7r~7r°7, 7r + 7r~7r + 7r~7, 
7r + 7r~7r + 7r°7, 7t + t/7) (|22ip in the mass range 0.6 < Mx d < 1.8 GeV, which covers 
approximately 50% of the total branching fraction. To reduce the uncertainty 
due to missing modes and phase space, X s j modes were also measured in the 
corresponding seven final states in the same mass range, where the first ir + was 
replaced with K + . The ratio of the two inclusive branching fractions is 

-J — — — ^r = 0.033 ± 0.013 ± 0.009, (54) 

B{B ->• X s 7j 

which is converted to |Vtd/Vt s | = 0.177 ± 0.043 ± 0.001, where the theory error 
does not include the effect due to the limited mass range. This result is also in 
good agreement with other determinations. 

The direct CP asymmetry for B — ¥ py can be as large as ~ —10% in the 
SM, whereas the time-dependent CP asymmetry is doubly suppressed due to the 
photon helicity and the cancellation of the CKM element Vtd- The latter appears 
in the mixing and in the b — > d penguin decay. However, the p° — > ir + ir~ decay 
provides clear vertex information for B° — > p 0r y. Both CP asymmetries have 
been measured by Belle (|219jl222p as 

A C p(B + -)■ p+7) = -0.11 ± 0.32 ± 0.09, 

A C p(B° -> p° 7 ) = -0.44 ± 0.49 ± 0.14, and (55) 

S C p(B° -t p°j) = -0.83 ± 0.65 ± 0.18. 

So far the results are consistent with null asymmetry. A nonzero direct CP 
violation may be measured earlier in B — > /ry than in B — > K*j. 

The isospin asymmetry (Eq. [35]) is also expected to be as large as ~ —5% 
inB-> pj. Belle measures A(pj) = -0.48 ±g;fg ±%^ (l2T9|) . and BaBar mea- 
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sures A(/ry) = — 0.43 JI022 -^ ^-10 (220); both measurements show a large isospin 
asymmetry. The average by HFAG is 

A(p7) = -0.4fi±{};& (56) 

A significant nonzero isospin asymmetry could indicate NP. 

5.5 Exclusive B — > K^£ + £r Branching Fraction 

Despite their small branching fractions, the exclusive decay channels B — > K^£ + £' 
have been measured efficiently with small background at Belle and BaBar, given 
that their final states are the same as those of B — > J/ipK^*' for which the B 
factories were designed. Here, K^ is one of K + , Kg, K* + and K*°, and £ + £~ 
is either e + e~ or /i + /i~. 

Electrons are identified by their energy deposit through an electromagnetic 
shower in the calorimeter. The minimum momentum is required to be greater 
than 0.4 GeV by Belle or 0.5 GeV by BaBar. The momentum of the bremsstrahlung 
photons that may be emitted by the electrons are added to their momenta if they 
are found near the electron direction. Muons have to reach and penetrate into 
the outer muon detectors and the minimum momentum is required to be 0.7 GeV 
by Belle or 1.0 GeV by BaBar. The dilepton mass regions around J/tp and ip(2S) 
are vetoed. 

The branching fractions, averaged over the lepton and kaon flavors and inte- 
grated over the dilepton masses, assuming the SM distribution over the vetoed 
J/ij) and ip(2S), were measured by Belle (I223I) and BaBar (|224p and have been 
averaged by HFAG as 

B(B^ K£ + £~) = (4.5 ±0.4) x 10~ 7 , . 

B(B -* K*£+£-) = (10.8 t\i) x lO" 7 . [ ' 

The results are consistent with SM expectations. At present, the irreducible 
form factor uncertainty in the SM calculations prevents these results from placing 
meaningful constraints on NP. 

A small subset of these combinations, B + — > K + fi + fi~ and B° — > K*° u. + u.~ , 
can be efficiently measured at hadron colliders. CDF has reported the most 
precise measurements of these modes (225). 

5.6 B — ¥ K^*H + £~ Asymmetries and Angular Distributions 

The direct CP and isospin asymmetries in B — > K^*>£ + £~ are also useful in 
the search for NP. The direct CP asymmetries are consistent with null values, 
A C p{B -► K+£+£~) = -0.01 ± 0.09 and A C p(B ->■ K*£+£~) = -0.07 ± 0.08 as 
averaged by HFAG. However, nonzero negative isospin asymmetries in the small 
<? 2 -region of B -)■ K£ + £~ (3.2a) and B -> K*£ + £~ (2.7a) have been reported by 
BaBar (3.9<r when combined). The corresponding isospin asymmetries by Belle 
are 1.4cr and 1.8c from zero, and are consistent with both BaBar 's results and null 
asymmetry. The isospin asymmetry combined for B — > K^*>£ + £~ and averaged 
by HFAG is 

Af M = -0.45 ±0.10. (58) 



The SM prediction is essentially zero at this level of statistics (see Section 14.2 
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Muon to electron ratios in B — > K^*>£ + £~ (Eq. [30]) are also measured by both 
Belle and BaBar. Results are consistent with the SM, and their naive averages 
are R K = 1.02 ± 0.18 and R K * = 0.88 ± 0.17. 

The four-body decay configuration of B — > K*£ + £~ — > Kir£ + £~ allows ex- 
traction of further information from the angular distributions of the final-state 
particles. The most interesting observables are the fraction of longitudinal polar- 
ization Fl from the kaon angular distribution (Eq. [4Tj) and the forward-backward 
asymmetry Afb from the lepton angular distribution (Eq. I42J) . Belle has mea- 
sured Fi and Afb in six bins of q 2 (223), whereas BaBar has done so in two 
bins (226). Current statistics are not enough to tell whether there is a zero- 
crossing point at low q 2 , although the results favor the case with no crossing, 
for which the sign of the Wilson coefficient Cj is flipped. Both results have pos- 
itive Afb for high q 2 (Fig. [6|), which sets nontrivial constraints on the Wilson 
coefficients. CDF has also measured Fp and Afb in the same six bins as Belle 
for B° — > K*°fi + fi~ events (|225p . The results are in agreement with Belle and 
BaBar. 

5.7 Inclusive B — > XJ + £r Branching Fraction 

The inclusive B — > X s £ + £~ branching fraction has been measured by Belle and 
BaBar using the sum-of-exclusive technique. The X s system includes final states 
with one kaon and up to four (two) pions that have masses up to 2.0 (1.8) GeV for 
the result by Belle (BaBar). Belle recently announced a preliminary result based 
on 657 million BB p27]) . and BaBar 's result is based on 89 million BB (|228j) . 
In Belle's new analysis, partial branching fractions are measured in bins of the 
X s mass, and then the total branching fraction is calculated as their sum. This 
method reduces the large systematic error observed in previous studies that arose 
from the strong X s mass dependence of the efficiency and the unknown fractions 
of exclusive channels B — > K^*'£ + £~. The measurement is still dominated by 
the statistical error and will be more precise in the future. Belle and BaBar 
reported branching fractions as B{B -> X s £ + £~) = (3.33 ±0.80 toil) x iO" 6 and 
B{B — > X s £ + £~) = (5.6±1.5±0.6±1.1) x 10~ 6 , respectively, which were averaged 
by HFAG as 

B(B -»■ X s £ + £-) = (3.66 ±° f 7 ) x 10^ 6 , (59) 

integrated over the entire subset of phase space with q 2 > 0.2 GeV, including the 
vetoed J ftp and ip{2S) regions. The results are in good agreement with the SM 
prediction. They strongly disfavor the case with the flipped sign of C-j (I179p . 

6 OUTLOOK 

Remarkably, the B factories have measured all the observables within the radia- 
tive and electroweak penguin decays at values that are consistent with the SM 
predictions. These measurements rule out 0(1) corrections to the SM and iden- 
tify the CKM theory as the dominant effect for flavor violation as well as for CP 
violation. The success of the simple CKM theory of CP violation was honored 
with the Nobel Prize in Physics in 2008. Theoretical tools and precision have 
significantly advanced during the past decade, and we are ready to challenge the 
SM if a clear deviation is found or to discriminate different NP scenarios if direct 
evidence is found at the LHC. 
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Also, the future offers great experimental opportunities in flavor physics. LHCb 
has finally started taking data and promises to overwhelm many B factory results, 
and ATLAS and CMS will also contribute to flavor physics. In the radiative and 
electroweak penguin decays, the most promising measurements are the angular 
analysis of B° — > ET*°/x + /i~ and the analysis of time-dependent CP asymmetry 
in B s — > 07; the latter measurement cannot be performed at the B factories due 
to the fast B s oscillation. However, the theoretically clean inclusive modes and 
many modes involving neutral particles like the 7r° can be pursued only at the 
e + e~ B factories. Two proposed super- B factories, Belle II at KEK and SuperB 
in Italy, would accumulate two-orders-of-magnitude-larger data samples. Such 
data would push experimental precision to its limit. 

Theoretical and experimental techniques are ready for such large data samples. 
The results provided by LHCb and the next-generation e + e~ B factories are 
eagerly awaited, as they may be the key to identifying physics beyond the SM. 
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Figure 1: Examples of radiative penguin decay diagrams (a) in the Standard 
Model and (b) beyond, (c) A penguin. 
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Figure 2: Renormalization-scale dependence of B{B — > X s ^f) in units 10 -4 at 
leading log (dotted lines), next-to- leading log (dashed lines) and next-to- next-to- 
leading log (solid lines). The plots describe the dependence on (left) the the low- 
energy scale u b and (right) the charm mass renormalization scale fj, c , from (j58[) . 
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Figure 3: Spectra in inclusive modes: (left) Cut in the photon-energy spectrum 
in B — > X s "f. (right) Differential B — > X s l + £~ branching fraction as a function 



of s 



q /m 



rn 



£+t 



-/mi, including the effect of charm resonances in the 



Kriiger-Sehgal method (solid line). For comparison, the dashed curve shows the 
same quantity obtained within a purely partonic calculation at next-to-leading- 
log precision, from Ref. (|113p . 
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Figure 4: Example of Mes (Mb c ) and AS for B° — > K*°~f by BaBar (from 
Ref. (HMD)- 
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Figure 5: Photon-energy spectrum for i? — )• X s 7, as measured by Belle (a) with- 
out lepton tag, (b) with a lepton tag, and (c) their average from Ref. (196). 
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Figure 6: Longitudinal polarization fraction, forward-backward asymmetry and 
isospin asymmetry of B — > K*l + l~ by Belle (right, from Ref. (223)) and BaBar 
(left, from Ref. ([226J)). The solid line shows the Standard Model predictions, and 
the other curves represent non-Standard Model extreme cases. 



